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ABSTRACT 

We report the discovery of TeV gamma-ray emission from the Type la supernova remnant (SNR) 
G120.1-I-1.4, known as Tycho's supernova remnant. Observations performed in the period 2008- 
2010 with the VERITAS ground-based gamma-ray observatory reveal weak emission coming from the 
direction of the remnant, compatible with a point source located at 00*^ 25™ 27. 0'', -1-64° 10' 50" 
(J2000). The TeV photon spectrum measured by VERITAS can be described with a power-law 
dN/dE = C{E/3A2 TcY)-'^ with F = 1.95 ± 0.51stat ± O.SO^ya and C = (1.55 ± OASstat ± OAlsys) x 
10~^^ cm~^s~^TeV^^. The integral flux above 1 TeV corresponds to ~ 0.9% percent of the steady 
Crab Nebula emission above the same energy, making it one of the weakest sources yet detected in 
TeV gamma rays. We present both leptonic and hadronic models which can describe the data. The 
lowest magnetic field allowed in these models is ~ SOfiG, which may be interpreted as evidence for 
magnetic field amplification. 

Subject headings: gamma rays: observations — ISM: individual objects (G120. 1+01.4, Tycho=VER 
J0025-f641) 



1. INTRODUCTION 
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The object G120.H-1.4, commonly called Tycho's su- 
pernova remnant (SNR), is the historical relic of a su- 
pernova which was first observed in 1572. Recent spec- 
tral a nalysi s of the light echo from the explosion ( Krause] 
et al. 2008 1 has confirmed the long-standing conjecture 



see, e.g., Ruiz-Lapuente 2004 ) that the event was a Type 



la supernova. 

At 438 years old, Tycho is young among Galactic SNRs 
and is well-studied at many wavelengths. It has a radio 
spectral index of 0.65 and a flux density at 1.4 GHz of 
40.5 Jy (Kothes et al. 2006| ). A hint of spectral curva- 



ture, consistent with nonlinear shock acceleration, and 
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a slight ly flatter radio spectrum ( 0.61) has also been re- 
ported ( Reynolds fc Enison]|1992[ ). Radio images show a 
clear shell-like morphology with enhan ced emission along 
the northeastern edge of the remnant ( Dickel et al.|1991 
Stroman fc Pohll[2009| ) 



-X.-ray images reveal strong non-thermal emission con- 
centrated in the SNR rim (Hwang ct al. 2002; Bamba 

al., 



et al. 2005 Warren et al. 



2005j |Katsuda et 



, , 2010| ). 

i'hin nlamentary X-ray structures in this region have 
been interpreted as evidence for electron acceleration 



Hwang et al. 2002 Bamba et al. 2005 Warren et al 



U05| . This IS supported by ra dio spectral tomography 



studies (Katz-Stone et al. 12000 ) and b y the detection of 



X-ray s up to energies ot 30 keV (e.g., iTamagawa et al 



2009 1 , which implies the presence of electrons up to at 
least ~ 10 TeV. 

The global expansion rat e of the remnant, measured 
at many wavelengt hs (Kam per fc van de n Bergh 1978 



Reynoso et al. 1997; Katsuda et al. 2010 ), is consistent 



with an object transitioning into the Sedov phase. How- 
ever, many of these studies have noted that the north- 
east quadrant of the remnant is expanding at a lower rate 
than the rest of the object. This has been attributed by 
some to interactions with a nearby high-density cloud, 
which has been studied through its H I and CO emissi ons 
( [Reynoso et"al]|1999| |Lee et al.|[2004l |Cai et al]|2009[ ). 

Distance estimates tor I'ycho have varied. Estimates 
which combine measurements of proper motion with 
shock velocities inferred fro m Hp. line widths tend to fa- 
vor distances near ~ 2.5 kpc (Chevalier et al. 1980 Smith 
On the other hand, the distance de- 



et ah 1991 1992) 



rived from the light echo spectrum ( Krause et al. 2008 1 
is 3.8^i['^ kpc, and HI absorption studies (Schwarz et al. 



1995 1 yield 4.5 ± 0.5 kpc. A recent X-ray study com 



bining ejecta v elocities with proper motion studies yields 

4.0 ± 1.0 kpc (iHayato et al. 120101). The apparent asso- 
THh. 



elation with a molecular CO cloud h as not produced an 
unam biguo us distance estimate (e.g.,|Lee et al.|2d04 Cai 



et al. 2009), due to the complex velocity held in the di 
rection ot the Perseus arm, where Tycho appears to be 
located. For the same reason, estimates based on the HI 
measurements ma y be considered controversial. Finally, 
Volk et al. (20081 derive a lower limit of 3.3 kpc, based 
on the previous non-detection of Tycho in TeV gamma 
rays. 

It is worth noting that many of the shorter distance 
estimates rely on assumptions about the post-shock pro- 
ton temperature's relationship to the sho ck velocity in 
the case of adiabatic expansion (see also Strom 1988). 
These assumptions may not be valid in the presence of 
ef ficient particle accel eration, as discussed, for instance. 



Helder et al.| ( 2009[ ). 

Claims ot sucn ethcient nuclear particle acceleration 



have in deed been made (e.g., Warren et al.|2005 though 
see also Lee et al.|20 10l , based on detailed studies of the 
shock front and contact discontinuity locations. Draw- 
ing on this result and on their own detailed simulations, 
Volk et al. ( 2008 1 claim that a detection of gamma rays 
from 'I'ycho would represent "incontrovertible" evidence 
of nuclear particle acceleration in SNRs. More generally, 
TeV gamma-ray observations probe the high-energy end 
of the underlying parent particle distributions and, espe- 
cially in conjunction with measurements at other wave- 



lengths, can help c onstrain the associa ted acceleration 
mechanics (see, e.g., Kaaret et al.||2009" ). 

Tycho has been observed several times at gamma-ray 
energies, with no credible detections previously reported. 
In the GeV regime, the source does not appear i n either 



the first Fermi LAT catalog (Abdo et al. 



3rd EGRET catalog ( [Hartman et al.||l999| Mj^per lim 



2010) or the 



its in the TeV range have been pres ented by the Whipple 
Collabo ration ([Buckleyet al.|1998| ), the HEGRA CoUab- 
oration (JAharonian ct ai. 2001') and the MAGIC Collab- 
oration (Carmoiia et ai. 20(J9j. This last limit is the most 
constraining, with a 3(f~(where a indicates standard devi- 
ations) flux upper limit of about 1.7% of the steady Crab 
Nebula flux above 1 TeV. All of the TeV limits assume a 
gamma-ray source located at the center of the remnant. 

In this Letter we report the discovery by VERITAS of 
TeV gamma-ray emission from Tycho's supernova rem- 
nant. 

2. VERITAS INSTRUMENT & OBSERVATIONS 

VERITAS is an array of four 12 m imaging atmospheric 
Cherenkov telescopes located at the Fred Lawrence 
Whipple Observatory in southern Arizona (1.3 km a.s.L, 
N 31°40', W 110° 57'). Each of the telescopes has a cam- 
era comprising 499 photomultiplier tubes covering a 3.5° 
total field of view. The array is sensitive to photon ener- 
gies between 100 GeV and 30 TeV, with an energy reso- 
lution of ~ 15% and an angular resolution per event of 
^ 0.1°. The array has been fully functional since Septem- 
ber 2007. In summer 2009, telescope Tl was relocated, 
impro ving the sensitivity and angular resolution of the 
array ( Perkins et al.|2009 ). Prior to this reconfiguration, 
VERIIAS was capable o: 



' detecting a point source with 
1% of the flux of the Crab Nebula at the 5a level in less 
than 50 hours. After moving Tl, such a source can be 
detected in less than 30 hours. 

VERITAS observations of Tycho spanned two epochs, 
from October 2008 to January 2009 and from Septem- 
ber 2009 to January 2010. All observations were taken 
in "wobble mode" , in which the telescopes are pointed 
to a position 0.5° offset from the source position. The 
offset direction is sequentially varied from run to run be- 
tween the four cardinal directions. Approximately equal 
amounts of data were acquired from each offset position. 
This technique allows for the collection of data with a 
simultaneous estimate of the background. From the pe- 
riod between October 2008 and January 2009, a total of 
21.9 hours of data were retained after quality cuts on 
weather conditions and operational stability. The mean 
zenith angle of these data is 35°. Between September 
2009 and January 2010, an additional 44.7 hours of data 
were analyzed after quality cuts. The mean zenith an- 
gle of this data set is 39° . Note that the second season's 
data set was taken with the new array configuration, and 
therefore benefits from the improved sensitivity. 

3. DATA ANALYSIS AND RESULTS 

The da ta were analyze d following standard procedures 
(see, e.g., Acciari et al.|2008) . The optimum gamma- ray 
selection criteria (cuts) depend upon the source flux and 
spectrum, which are not known a priori. We therefore 
define and apply two sets of cuts, appropriate for a rea- 
sonable range of source properties, and account for the 
statistical trials incurred in this process. The results pre- 
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Fig. 1. — VERITAS TeV gamma-ray count map of the region 
around Tycho's SNR. The color scale indicates the number of ex- 
cess gamma-ray events from a region, using a squared integration 
radius of 0.01 deg^ for the 2009/2010 data and 0.015 deg^ for the 
2008/2009 data. The centroid of the emission is indicated with a 
thick black cross. Overlaid on the image a re X-ray conto urs from a 
Chandra ACIS exposure (thin black lines; [Hwang et"aT.|2002j ) and 
^^CO emission (J=l-0) from the high-resolution FCRAO Survey 
(magenta lines; [Heyer et al.|19 98|). The CO velocity selection is dis- 
cussed in the text. Ihe VJiKiTAS count map has been smoothed 
with Gaussian kernel of size 0.06°. The point-spread function of 
the instrument (see text) is indicated by the white circle. 

sented here required that at least 3 of the telescopes in 
the array recorded an image with more than 800 digital 
counts (~ 150 photo-electrons, corresponding to an en- 
ergy threshold of 800 GeV) . These images were then used 
to reconstruct the air shower properties. Gamma-ray 
events were selected through cuts on the mean reduced 
scaled width and mean red uced scaled length parameters 

which were required to fall 



(Krawczynski et al. 



2006) 



between -1.2 and 0.5, and on the square of the angu- 
lar distance from the test position to the reconstructed 
arrival direction of the shower, 6"^. For the 2008-2009 
data set, we required 9^ < 0.015 deg^, and for the 2009- 
2010 data set, we required 9^ < 0.01 deg^ because of the 
improved angular resolution of the array. The ring back- 
grou nd model was used to estimate the background (see, 
e.g., Berge et aL]|2007[ for a description). 

This analysis jaroduced an excess which is significant 
at the 5.8(7 level (pre-trials). This is the peak excess 
found in a blind search region with sides of length 0.26° 
- roughly twice the diameter of the radio remnant. A con- 
servative a priori trials factor w as determined by tilin g 
this area with square 0.04° bins (Aharonian et al.||2006), 
and additionally accounting for the two sets o± applied 
cuts. This results in a 5. Oct post-trials significance. 



3.1. Morphology 

The morphology of the source was investigated by bin- 
ning the uncorrelated acceptance-corrected map of excess 
event counts. Bins of size 0.05° were used to provide suffi- 
cient statistics for fits to source models. The map is com- 
patible (x^ = 508; ndf = 438; Prob = 1.2%) with a point 
source located at 00^ 25" 27.0", -f64° 10' 50" (J2000) 
and hence we designate the object VER J0025-I-641. This 
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Fig. 2. — Differential gamma-ray photon spectrum of Tycho 
as measured by VERITAS. The error bars represent la statistical 
errors only. The solid red line shows the results of a power-law fit 
to the VERITAS data. The lower panel shows the residuals of the 
data fror n this fit. T he dashed blue line represents the hadronic 
model of I Volk et a!] pOOSi for gamma-ray emission from Tycho, 
scaled to 4 kpc. i ne details of the model and of the analysis are 
discussed in the text. 

position is derived from a simple symmetric Gaussian 
fit with a width fixed at the instrument point-spread 
function (i.e., the 68% containment radius for photons, 
0gg%= 0.11°). While other source functions (e.g., an 
offset asymmetric Gaussian) may provide a marginally 
better fit (perhaps hinting at a more complex underlying 
source morphology) a likelihood ratio test shows that the 
extra degrees of freedom are not statistically justified in 
this data set. 

As shown in Figure 1, the center of the fit position is 
offset by 0.04° from the center of the remnant. The sta- 
tistical uncertainty in this location is 0.023°, while the 
systematic uncertainty resulting from telescope pointing 
accuracy is 0.014°. We note that the derived centroid po- 
sition depends also on the source shape assumed for the 
fit. Future observations will allow more detailed study of 
the source morphology. 

3.2. Spectrum 

The differential photon spectrum between 1 and 
10 TeV is shown in Figure [2j This spectrum is gen- 
erated from the complete data set after quality selec- 
tion. The shape is consistent with a power law dN/ dE = 
C(S/3.42 TeV)-r with Y 
C= (1.55 ±0.43,tat± 0.47 



1.95±0.51,tat±0.30sy. and 



' sys 



X 10- 



-14 —2 

cm s 



TeV" 



where the systematic error on the flux is dominated by 
uncertainty in the energy scale. The of the flt is 0.6 
for 1 degree of freedom. The integrated flux above 1 TeV 
is (1.87±0.51stat) X 10-13 cni-2 g-i, about 0.9% that of 
the steady Crab Nebula flux above the same energy. 

4. DISCUSSION 

Figure [T] shows the TeV gamma-ray image of Tycho's 
SNR. The color scale indicates the number of excess 
gamma-ray events in a region, using a squared integra- 
tion radius of 0.015 deg^ for the 2008/2009 data and 
0.01 deg2 for the 2009/2010 data. The map has been 
smoothed with a Gaussian kernel of radius 0.06°. Over- 
laid on the image are X-ray contours fr om a Chan dra 
ACIS exposure (thin black lines; , Warren et al.||2005 ). A 
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Fig. 3. — Radio, (nonthermal) X-ray, and VHE gamma-ray emis- 
sion from Tycho's SNR, along with models for the emission. The 
upper panel shows a lepton-dominated model while the lower panel 
shows a model dominated by hadrons. In each, the IC emission 
corresponds to the (cyan) long-dashed curve while the pion-decay 
emission corresponds to the (magenta) short-dashed curve. The 
solid curve at high energies is the sum of these components; at 
lower energies it corresponds to the synchrotron emission. See text 
for discussion. 

contour map of 115 GHz line emission associated with the 
molecular ^^CO (J=l-0) transition, from the 14 m tele- 
scope of the Five College Radio Astr onomy Ob servatory 
(FCRAO), is shown in magenta (Heyer et al.|19 98; Taylor 
et al.||2003| . Following the analysis ot |Lee et al. (2004), 
we have integrated over the velocity range -68 km s^^ to 



ence from the molecular cloud. Two versions are consid- 
ered - one in which the TeV emission is dominated by 
leptonic (inverse Compton (IC) scattering; upper panel) 
processes, and one in which it is dom inated by hadroni c 
(pion decay; lower panel) processes (Slane et al. 2010). 
Here we have assumed particle spectra ot the torm: 



dN 
dE 



(1) 



with the spectral index a being fixed to the same value 
for the electrons and protons, but with the cutoff energy 
Ec being allowed to differ, as expected for loss-dominated 
distributions. We note the crucial point that, in both the 
leptonic and hadronic models, the total relativistic par- 
ticle energy is dominated by hadrons, and it represents a 
significant fraction of the total su pernova kinetic e nergy, 
estimated to be 1.2 x 10^^ ergs ( Badenes et al.| [2006). 
Indeed, in both cases, the energy density ot hadronic 
cosmic rays is likely sufficient to modify the su pernova 
shock dynam ics, supporting the conclusions of IWarren 
et al. ( 2005 ) . The radio spectrum requires a 2.2, and 
we have assumed an ambient density no = 0.2 c m~^. 



based on u pper limits from X-ray measurements (Kat- 
suda et al.||2010) . The distance is set to 4 kpc. 

I'he ramo data in Figure |3] are compiled from Re ynolds] 
& EUison (1992). The X-ray data (shown in blue) rep- 
resent the unfolded spectrum between ~ 4 and 10 keV 
extracted from a deep Suzaku observation of Tycho. Fea- 
tures from Fe-K and weaker line emission have been re- 
moved. The TeV data points are from the VERITAS 
results reported here. 

For the lepton-dominated scenario, we find the nor- 
malization and magnetic field required to reproduce the 
radio and X-ray emission by synchrotron radiation. The 
gamma-ray emission is then generated by IC scattering 
of cosmic microwave background photons. (The impact 
of adding an additional - even maximal - cont ribution 
from the dust-IR emission ( Douvion et al.pOOl ) was in- 
vestigated and found to make only a small difference in 
the overall IC spectrum.) The associated magnetic field 
is ^ 80/iG, with a ^ 15% uncertainty. Assuming an 



electron-to-proton number ratio k, 



ep 



10- 



as observed 



in the local cosmic-ray population, the total particle en- 



50 km s-\ revealing a cloud possibly interacting with ergy (dominated by protons) is 1.8 x 10^° erg, and the as 



the northeast quadrant of the remnant 

As can be seen from Figure [l] the peak of the gamma- 
ray emission is displaced somewhat to the northeast of 
the center of the remnant, towards the CO cloud. While 
this is provocative in the contex t of hadronically-induce d 



sociated gamma-ray emission from tt^ decay is negligible 
in the TeV band, as shown in the upper panel of Figure[3j 
The magnetic field required for this model is somewhat 
lower than the conventionally accepted ~ 200 — SOO^tG 
value that many derive (e.g. 



gamma-ray emission (see, e.g 
the statistical significance of the displacement is weak 
Nevertheless, this is the general morphology expected in 
the case of a shock/cloud interaction leading to gamma- 
ray emission, which is se en in several other remnants (see, 



Aharonian et al.||1994|), 2008 Marcowith fc Cassell2010 1, tho ugh it is within the 



range found by Warren et aE 12005) 



Ballet 2006 Volk et al 



Hewitt et al. 2009). On the other hand, OH maser 



e.g., ^ ^ 

emission, a telltal e sign of such interactions (Wardle & 
Yusef-Zadeh 2002) has not been detected from this rem- 
nant (Frail et al. 1996). A catalog search within our error 
box reveals no likely gamma-ray emission candidates at 
other wavelengths, so we tentatively associate the source 
with Tycho. 

In Figure [3] we show simple model fits to the broad- 
band spectrum of Tycho, generated assuming no influ- 



For the hadron-dominated scenario, we adjust the nor- 
malization of the proton spectrum to reproduce the ob- 
served TeV emission through Tr^-decay (Figure [sj lower 
panel). The normalization of the electron spectrum is 
then reduced to a level at which the hadron-induced 
gamma-ray emission dominates, yielding a Kgp ~ 4 x 
10~*. The associated magnetic field required to repro- 
duce the synchrotron emission is ~ 230fj,G. This value 
can be considered a lower limit, since a smaller Kep will 
require a larger magnetic field. 

For this model, we find a total particle energy of 
~ 8 X 10^° erg. While this is perhaps uncomfortably 
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large, it is based on emission involving only the mean 
gas density of the remnant. Any possible contribution of 
target material from the cloud would reduce the energy 
requirement. This could also be achieved if we loosen the 
assumption of a single power-law particle spectrum. 

Our hadronic m odel is broadly consistent with that of 
Volk et al. (2008), who employ a nonlinear kinetic par- 
ticle acceleration model to derive a hadronically-induced 
flux as a function of distance. As shown in Figure [2 
with a source distance of 4 kpc, an ambient density of 
~ 0.2 cm~'^, and a magnetic field of ~ 350/iG, their 
model provides a reasonable fit to our data. However, by 
scaling the flux to improve the fit, a distance estimate of 
3.8 kpc can be obtained. 

Finally, it is worth noting that the lowest magnetic field 
allowable in either of our models is ~ 80/iG. This value 
is not only well above the typical ^ 3/iG fields of the 
interstellar medium, it is significantly higher than expec 
tations from shock-compression of that field in moderate 
densi ty environments (see 



e-g-, 



Ellison & Vladimirov 



2008|. Hence, our detection of gamma rays from I'y- 
cho's SNR represents additional evidence for magnetic 
field amplification in this remnant. While this conclu- 
sion is dependent on the validity of our one-zone emission 
approximation, it is not subject to large systematic un- 
certainties arising from the choice of electron spectrum, 
since, for a ~ 80/zG field, the Suzaku and VERITAS data 
probe the same underlying electron energies. 

5. CONCLUSION 

VERITAS has detected weak TeV gamma-ray emis- 
sion from the SNR G120.1-I-1.4, also known as Tycho's 
supernova remnant. The total fiux from the remnant 
above 1 TeV is 0.9% that of the Crab Nebula, mak- 
ing it one of the weakest sources yet detected in TeV 
gamma rays, and only the second confirmed Type la SNR 
gamma-ray emitter. SN1006, another Type la remnant, 
was observed by H.E.S.S. to have two distinct regions 
of weak 1% Crab) gamma-ray emission, one at the 
northeastern edge of the ~ 0.5° diameter remnant and 



one at the southwestern edge ( Acero et al. 2010 ). Both of 
these regions are highly correlated with non-thermal X- 
ray emission and have spectra that are compatible with 
power-laws of index F ^ 2.3. For SN1006, neither a lep- 
tonic nor a hadronic origin for gamma-ray emission can 
be firmly eliminated. 

The photon spectrum of Tycho can be described by a 
power-law with differential index 1.95 ib O.Slst at =1= 0.30s^s 
and is compatible with the hadronic model of |Volk et al.| 
(20081, scaled for a distance of 3.8 kpc. We present an- 
other hadronic model which can also reproduce the data, 
although it requires a large energy content of cosmic rays. 
A leptonic external IC model also provides a tolerable fit 
to the data, though it requires a magnetic field somewhat 
lower than generally is accepted for Tycho. Notably, the 
lowest magnetic field allowed in either model is ~ 80/xG, 
which may be interpreted as evidence for magnetic field 
amplification. The morphology of the emission is com- 
patible with a point source, with the peak of the emission 
possibly offset from the center of the remnant. Detailed 
spectral and spatial studies will be possible with a deeper 
exposure. 
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